Charged Plate Beyond Mean-field: One-loop Corrections by Salt-density Fluctuation 

William Kun£| 

Department of Physics, Syracuse University, Syracuse, New York 13244 



O 
O 

(N 



oo 



o 



I 

C 

o 
o 



> 

m 
O 

O 



T3 

O 
o 



13 



A.W.C. Lau 

Department of Physics, Florida Atlantic University, Boca Raton, Florida 33431 

(Dated: February 6, 2008) 

We present an exact field-theoretic formulation for a fluctuating, generally asymmetric, salt den- 
sity in the presence of a charged plate. The non-linear Poisson-Boltzmann equation is obtained as 
the saddle-point of our field theory action. Focusing on the case of symmetric salt, we systematically 
compute first-order corrections arising from electrolytes fluctuation to the free energy density, which 
can be explicitly obtained in closed form. We find that for systems with low to moderate salt den- 
sity, fluctuation corrections to the free-energy depends sensitively on the salt concentration as well 
as their charge valency. Further, we find that electrolyte fluctuation leads to a reduced electrostatic 
repulsion between two point-charges when they are close to the charged plate. 

PACS numbers: 82.70.-y, 61.20.Qg, 05.70.-a 



I. INTRODUCTION 

Across systems from astronomical scale to cellular 
level, electrostatic interaction plays an important role 
in their characterization by thermal properties. For 
these systems, electrostatic interactions are modulated 
by charged ions interacting via Coulomb and hard- 
core potentials. Examples include high-density plasma 
found in stellar interior and solar photosphere Q, 0, 
charged colloids found in such industrial applications 
as coating materials and ceramic precursors ]1 0] , and 
the self-assembly of charged surfactants into lamellar 
phases d |(| . Water itself can be considered a weak elec- 
trolytic solution at room temperature, and as a result, 
considerations of electrostatic interactions figure promi- 
nently into many biological problems involving macroions 
in aqueous solutions |3- As examples, the macroions may 
be charged cell membranes or stiff polyelectrolytes such 
as DNA and polymers d ||- Often, the interplay be- 
tween thermal fluctuation and electrostatic interaction 
in these charged systems determine their chemical affin- 
ity, structural and functional properties. To ensure their 
stability, there must also be overall charge neutrality. In 
experimental conditions, the charged particles in aque- 
ous solutions often exist in the form of dissociated cations 
and anions. Due to the long-range nature of Coulomb in- 
teractions, these charged particles can become correlated 
and drive the vastly complex behaviors observed in these 
systems. 

Quantitatively, it is assumed that the cationic and an- 
ionic densities in these charged systems are conserved and 
obey the Boltzmann distribution in equilibrium [ij. The 
conventional description then relies upon the mean-field 
approach in which one averages out the fluctuation in 
the distribution of salt ions. Combining Gauss' law and 
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the Boltzmann distribution of the salt ions then leads to 
the well-known Poisson-Boltzmann equation [lflj . The 
Poisson-Boltzmann equation provides a good quantita- 
tive description of most experimental systems and simu- 
lations [3, uM smce m general each charged ion simulta- 
neously interacts with many other ions in solution; one 
can reproduce the interaction as an effective mean-field 
experienced by each ion throughout the system. Despite 
the straightforward derivation of the Poisson-Boltzmann 
equation, exact solutions are only known in a few cases, 
when one imposes further simplifications such as assum- 
ing only the presence of counter-ions, restricting to prob- 
lems possessing planar 0, 0, 0] or cylindrical geome- 
tries flCI Il5j | , and linearization of the Poisson-Boltzmann 
equation in the appropriate regime [Hif . 

However, it is well known by now that the mean- 
field approach of the Poisson-Bolzmann equation works 
only in the low-charge regime 

HE III El El- There 
has been great interest to go beyond the mean-field de- 
scription of the electrostatic systems, since the mean- 
field nature of the Poisson-Boltzmann equation neces- 
sarily precludes correlation effects. The need to con- 
sider correlation arises, for example, in such problems as 
DNA condensation with multivalent cations 0, 0, 0] 
where the macroions are highly charged and the cor- 
relation effects between ions become important. Phe- 
nomena such as char ge i nversion and like-charge attrac- 
tion [53, HI El S Hi IH have been reported in 
these cases. Much experimental effort has been spent to 
confirm these predictions. Experimental techniques in- 
volving manipulating optical tweezers in confined geome- 
tries 1271 l28l I29I. interaction measurements in equilib- 
rium |30ll3lLl32ll33ll3^ | . and digital video microscopy [3{| 
have been used to measure attraction between like- 
charges in various systems. 

Aside from the importance of correlation effects which 
the mean-field approach, by construction, fails to incor- 
porate, it has also been shown theoretically that within 
the framework of the Poisson-Boltzmann theory, the in- 
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teraction of two identical charged particles is always re- 
pulsive [3(| EtJ. Modifications have been made to the 
Poisson-Boltzmann framework to incorporate the effects 
of fluctuation and correlation HI ES EE El El And 
attraction between highly-charged macroions and coun- 
terions have been demonstrated via Monte Carlo simula- 
tions |43| . numericals methods of hypernetted chain in- 
tegral equation [44. |45j and the local-density-functional 
theory |46l |47| . Netz and Orland have formulated the 
problem of fluctuating counterion distribution in the 
presence of a fixed charge distribution in the language 
of field theory [U EE EE EE HH1 • In particular, they have 
computed corrections to the electrostatic potential, coun- 
terion density, and interaction between charged particles 
in the presence of a fixed charge distribution in both the 
Poisson-Bolzmann limit and the strong coupling limit. 
They have also applied field-theoretic methods to obtain 
the first few coefficients of a systematic low-energy ex- 
pansion of the free energy for these charged systems as 
well as to investigate issues regarding charge renormal- 
ization and fluctuation-induced modification of the static 
van der Waals interactions of salt ions 0, EE EE EE Ell • 
The problem of a charged plate in salt solution has previ- 
ously been looked at by Dean and Horgan in the context 
of soap films 

El EE El E3 They too have formulated 
the problem in a field-theoretic language. However, their 
results have been obtained only after approximating the 
action to quadratic order in all terms. And their calcu- 
lations involve mainly the total pressure exerted on the 
film which does not require an explicit computation of 
the free energy of the system. 

In this paper, we extend previous works in two impor- 
tant ways. We systematically compute one-loop correc- 
tions, without further approximation, to the interaction 
between two charged ions near the plate due to their 
fluctuation. We also obtain an explicit expression for the 
one-loop free energy. Similar in approach to the previ- 
ous works, we formulate our problem in a field-theoretic 
language that naturally incorporates the effect of fluctu- 
ation and provides a recipe for the direct computation of 
higher-loop corrections to various quantities. We show 
that the nonlinear Poisson-Boltzmann equation can be 
derived as the saddle-point of the full field-theory action. 

In what follows, we will first establish the field- 
theoretic formalism and derive the Poisson-Boltzmann 
equation from first principles. The main steps involve 
writing down the microscopic Hamiltonian for the single 
charged plate in a generally asymmetric salt solution and 
rewrite the resulting partition function in terms of field- 
theoretic variables. This will be done in Section II of this 
paper. We will then specialize to the case of symmetric 
salt ions and briefly summarize results from the mean- 
field approximation and discuss the relevant physics for 
these salty systems in Section III. In the subsequent sec- 
tion, we will systematically consider corrections to the 
mean-field predictions when we incorporate the effect of 
fluctuating salt-ions densities . We will first expand the 
action around the saddle point and, using techniques of 



Green's function and functional determinant, compute 
one-loop correction to the free energy density. We will 
also consider the interaction between ions in the presence 
of the charged plate and their self energies in section IV. 
Finally, we summarize and conclude in Section V. 



II. FIELD-THEORY FORMULATION 

In this section, we formulate the problem of a charged 
plate in salt solution in a field-theoretic language and es- 
tablish notations. The field-theoretic approach has been 
previously established in the study of s yste ms contai ning 
counterions and symmetric salts \U. \48L l49l l50l l5ll 1521 153 . 
l54l E3 • We hereby extend the formulation to the general 
case of assymmetric salts. We consider a system that con- 
sists of N_ point-like particles of charge — Z_ |e| and N + 
point-like particles of charge +Z + \e\, as well as a single 
charged surface of constant charge density. Without loss 
of generality, we assume the surface to possess a negative 
charge density cr(x) = — |e| n 5(z) — — |e| n/(x); results 
for a positively charged surface can be obtained with min- 
imal modification. The charge neutrality condition can 
then be written as 



ZeN_ + ZeN: 



\e\ n Q A 







(1) 



where A is the area of the plane. Strictly speaking, the 
overall neutrality of the system must imply the relation of 
N + > A_ for the case of negatively charged surface; both 
the counterions and salt cations contribute to the overall 
distribution of positive charges. In what follows, we will 
not make the distinction between the counterions and 
cations; we will simply refer to both species as counter- 
ions while the salt anions as coions. The energy of the 
system can be simply written as: 



(3E N - Z 2 _l B + Z 2 + l B - 

j>k l x j < I j>k I x j < 

N + JV_ 



N- 



j=lfc=l \ X j X fc 



+z_X>(x7)-z + x;#xt) 

3 3 

where we have defined 



(2) 



= h 



j ™/( x ') 



(3) 



as the "external" field induced by the presence of the 
charged plate, and the Bjerrum length Ib = e 2 /ekT de- 
scribes the length scale where electrostatic interaction 
dominates thermal fluctuation. At room temperature, 
Ib ~ 7A for water. It is natural to define a number 
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density operator ^(x 1 * 1 ) = X)i=i <5( x± — x i) for both the associated with the interaction energy in Eq. J5J as fol- 
counterions and coions. In terms of these density oper- lows: 
ators, we can construct the canonical partition function 



1 N+ r r J + 1 i N - r r a - 1 N - N + 

n [/ if io n [/ «p ^-no) + ^iv + y (0 ) - z_ e *w + ^+ e *w 



dxdx' |Z +/ 3(x + ) - Z_p(x")J V(x,x') Z+p(x! + ) - Z-p(x! ) 



(4) 



where V(x, x') = | x _ x /| is the Coulomb interaction be- 
tween ions, and as such V(0) represents the Coulomb 
self-energy. In the form written in Eq. (0} , the Coulomb 
interaction acts on the effective charge density p c — 
Z + /5(x + ) — Z_/5(x~) which is physically modulated by 
the spacial distribution of each ion species and their re- 
spective valence. From now on, we will assume that the 



total number of ions (counter-ions and salt) in solution 
is conserved, i.e. N + + N- = N. 

To elevate Eq. Q to a partition function of field vari- 
ables, we borrow techniques established in field theory 
and make use of the following resolution of the identity 
operator 



i = Jv P ± 5(p ± -p(x ± ))= J Vp ± Vil> ± expji J dxV ± (x) [^(x) -^(x*)] 



(5) 



Upon insertion of Eq. (JSJ) into Eq. @ and simplifying, terms of field variables p ± (x) and V ;± ( x ) : 
we can now rewrite our canonical partition function in 



= J Vp+Vp-V^V^- exp j-i J dxdx' [Z+p + (x) - Z-p~(x)] F(x,x') [Z+p+(x') - Z_p"(x')] 
+i J dx[^+(x)p+(x)+^-(x)p-(x)]J 



dx 4 



= Z;fV(O) + Z + 0(x+)- l </>+(x+) 



1 

nJ. 



dx 



,Z?_V(0)-Z-.4(x-)-iil>- (x~) 



N- 



(6) 



In the form written in Eq. © , the field-theoretic version 
of the canonical partition function Zn + ,n_ [p ± , ip^ , cp] is 
rather cumbersome. To bring the expression to a more 
manageable form, we will introduce a chemical poten- 
tial (in unit of UbT) and consider the grand canonical 
ensemble: 

oo 

2^-= E (^ + f + (e'-f- Z N+ , N _ . (7) 

JV_,JV + =0 



We also explicitly perform the functional integral over 
field variables p ± (x) in Eq. ©. The operation involves 
performing a functional Gaussian integration followed by 
another integration over a delta functional which imposes 
the condition that = — t|— . Defining a new field 

^ = 5 ( z %~ )' we can explicitly write the grand- 
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canonical partition function in Eq. J7J) as follows: where Ao is the usual normalization constant and 



7V Jvfe- S ^^ 



(8) 



S[% 



j*(x) (-V 2 ) *(x) - 



4ttZ 



B Z^(0)+Z + «x)+tZ + *(x)+ /J+ _ ^ nl B z 2 _ V(0)-Z _^(x)-tZ -*QQ+ 



(9) 



where we have substituted the explicit expression for — V 2 (5(x — x'). We can then bring the grand-canonical 
the inverse of the Coulomb operator 4-7rF~ 1 (x, x') = partition function in Eq. @ into the standard form 



*(x) (-V 2 ) *(x) 



^+ Z+[i*(x)+*(x)] _ 

2 2 



Z_[i*(x)+0(x)] 



(10) 



by defining the quantity Ib, which has the dimension of 
length and is simply a rescaling of the Bjerrum length; 
the quantities n± which will play the role of inverse 
screening length of the counterions and coions, respec- 
tively; and the quantities which can be interpreted as 
the bulk densities of the counterions and coions, respec- 
tively, in solution: 



is 

,.2 



Pa 



A-kI b , 

e Z 2 ± V(0)+^ ± 
^3 



(11) 

(12) 
(13) 



Let us define </>±(x) = Z±<fi(pc). Based on our construc- 
tion of the partition function, the average counterion 
and coion densities (p± (x)) can be readily related to the 
grand-canonical partition function: 



<P±(x)) = 



SlogZ. 



<50±(x) 



(14) 



Extension to higher-order functional derivatives of the 
partition function with respect to the external field </>(x) 
brings about higher-order cumulant correlation functions 
in the standard way. From this point on, we will limit 
ourselves to consider only the case where both the coun- 
terions and coions have the same valency Z± = Z and 
the same chemical potential fj,± = fi. We note that 
the chemical potential in the grand canonical poten- 
tial is not uniquely defined, since the transformation 
(x) — > (x) + i fig and ju — > fi — po leaves the free 
energy invariant. We also note that the "external" field 



(x), determined by 
-V 2 c 



(x) = 47r/ B n/(x) , 



(15) 



contains an arbitrary constant. To make connection with 
the standard Poisson-Boltzmann equation, we determine 
the saddle point for the action in Eq. IjlOfl. with the 
conditions that k± = kq and Z± = Z, via the Euler- 
Lagrange equation: 



V 2 ^*o(x)] 



Zk e Z[j* o (x)+0(x)] 



Zk e -Z[i*o(x)+0(x)] 



(16) 



To bring Eq. (|16|l into a more familiar form, we define a 
shift in the field variable: <p(x) = Z [i^o(x) + <^>(x)]. In 
terms of the new field ip, Eq. (|16[) becomes 

V 2 [<^(x)] - k 2 sinh[<p(x)] 



-4irZl B n f (yL) , (17) 

where Z 2 . Eq. (|17|l is the Poisson-Boltzmann 

equation for a charged plate in salt solution that contains 
both counterions and coions. Our reinterpretation of the 
Poisson-Boltzmann equation as a saddle-point of an un- 
derlying field-theoretic action affords us a systematic way 
to include fluctuation effects and go beyond mean-field 
theory. 



III. MEAN-FIELD THEORY 

As stated, closed-form solutions exist for the Poisson- 
Boltzmann equation only for the planar geometry [T3 . IL3L 
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FIG. 1: A single charged plate in salt solution. The system is 
symmetric with respect to z — 0. Without loss of generality, 
the charged plate is assumed to have a positive surface-charge 
density a, with counterions (solid circles) and coions (open 
circles) present on both sides of the plate. 



Il4j and cylindrical geometry 0, 0| . In cases of other 
geometries, one has to resort to further approximations. 
For the case of spheres, one can apply the Derjaguin 
approximation when the interaction potentials between 
spheres are much larger than the range of the poten- 
tial |Toj . Under this condition one may simply neglect 
curvature effects at the lowest order and adapt results 
from the case of fiat plates. Thus, there is an intrinsic 
usefulness in understanding the solution of the Poisson- 
Boltzmann equation for even the simple case of planar 
geometry. 

The solution of the full nonlinear Poisson-Boltzmann 
equation for a single charged plate in salt solution is well 
known in literature. We will provide here a brief sum- 
mary of results and the relevant physics of the system. 
First, we note that the single charged-plate system is a 
one-dimensional problem whose solution for < z < oo 
is related to solution for — oo < z < by symmetry (Fig. 
^) . To solve Eq. I|17|) with the delta-function source term, 
we simply consider the associated homogeneous equation 
and match the following boundary conditions imposed by 
the charged surface: 



dtp (z) 



dz 
dip (z) 



dz 





Z(T S t_B 

2e 



2 
A 



(18) 
(19) 



where A = 



is the Gouy-Chapman length. The 



boundary condition at the plate in Eq. i|19l) is found by 
integrating the Poisson-Boltzmann equation (|17|l over the 



interval — e < z < e and taking the limit as e — > 0. We 
have also introduced a useful length scale A in Eq. H19|) 
for future convenience. The solution for the mean-field 
electrostatic potential (fo(z) can be readily found: 



ipo(z) = -4 tanlT 1 e -"(W+*o) 



(20) 



where the constant of integration zq is completely deter- 
mined by the relation 



k\ = sinh kzq 



(21) 



For this problem, we can construct two dimensionless pa- 
rameters that completely characterize the system. The 
first useful dimensionless quantity encodes informa- 
tion about the charge density of our plate. Given a highly 
charged plate in salt solution, ^ 3> 1. The second di- 
mensionless quantity kA, on the other hand, yields in- 
formation on the charge densities of the salt ions. From 
Eq. H20(l . we see that k _1 represents the screening length 
of the ions. Since a highly charged background of ions 
implies stronger screening, the condition that kX 3> 1 in- 
dicates either a dense distribution of moderately charged 
salt ions or high valency from the ions themselves. Con- 
versely the limit of nX <C 1 corresponds to the case when 
only counterions are present. In what follows, we will see 
that these two dimensionless combinations of physical pa- 
rameters also provide a convenient choice when express- 
ing some of our more complicated results, in a way that 
would highlight their underlying physics. Meanwhile us- 
ing the mean-field electrostatic potential po given by Eq. 
(|20ll . we can compute the corresponding charge density 
for the ions as well as the free-energy density at the mean- 
field level. First, the counterion and coion charge densi- 
ties, as a function of distance z from the charged plate, 
can now be written down in a straightforward manner: 



Pa 



I ± e -K(z + Z ) 
1 ip e -K(z + Z ) 



(22) 
(23) 



Our definition for the salt-ions charge densities in Eq. 
(|22|l is consistent with the form of fo{ z ) with the con- 
vention of a plus sign. The Hclmholtz free-energy density 
is given by 

0F = J dxp_(x) [logp_(x)a 3 -l] 

+ / dxp+(x) [log p+(x) a 3 - l] 

+ -i-|dx[V^(x)] 2 , (24) 



where a is the typical ionic size of the charges. The first 
two terms are the entropy of the counter- and co-ions, 
while the last term is the electrostatic energy of the sys- 
tem. The free energy (3F can in general be divided into 
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a volume contribution fy and an area contribution f a'- 
(3F = f v V + f A A. Upon substitution of Eqs. and 
(I23JI into Eq. (|24|l and further evaluation, we can find the 



Pfv 



The plot of Eq. (ffi))) is shown in Fig. J5J. For simplicity, 
we assume that a 3 ~ Ib\ and that the free-energy den- 
sity is units of £bX. In the limit that n — > 0, the total 
free energy of the system based on Eqs. I|25l) and H2tj|) 
correctly reduce to that for a system with counterions 
only, which has the form of the free energy of an ideal 
gas. 
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4 




for the case of symmetric salt ions and A^(x) = \&(x) — 
\l/o(x). By definition, the saddle point satisfies the mean- 
field equation, and the linear term in Eq. I|27(l van- 



mean-field free-energy densities fy and f a as a function 
of x = k\: 



(25) 



(26) 



FIG. 2: Plot of the effective mean-field free-energy density as 
a function of x = kA. The free-energy density is measured in 
units of £bX. We assume that a 3 ~ £bX 2 . 



IV. FLUCTUATION EFFECTS 



With the field-theoretic formulation in place for a sin- 
gle charged surface in salt solution, we can now systemat- 
ically investigate effects of fluctuating salt-ion densities. 
In what follows, we will first describe the general method- 
ology of our calculations before presenting details of our 
results in following subsections. To go beyond mean-field, 
we first expand the action in Eq. © about the saddle 
point 



• ,(27) 



(28) 



ishes accordingly. Since the quadratic term is Gaus- 
sian in nature, we can formally write the partition func- 
tion of this quadratic action as a functional determi- 



Po (log po a 3 - 1) 
2 



e B x (x + V?+T) i B x 



log 



2 
Jb~X 



- \A + x 2 ^j log { 



2a 3 



-T-T-v/l + £ 2 log ( 7= + 1 

^bA Wx 2 + 1 + 1 



5[*,0] - 5[* o ,0]+ / d 3 



5H 



<5*(x) 

where TL is the Hamiltonian given by 

1 fl 



A*( x ) + i / d 3 x / d 3 y d ' H 



A#(x)A#(y) 



H = -f(x)(-V 2 )*(x)-K 2 cosh[Z[^(x) + (/ ) (x)]] 
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nant [56(. To compute this determinant, we define the 

kernel K(x,y) = g^My) ' and we observe tnat its °P~ 
erator inverse is just the Green's function G(x, y) defined 
by: 



J dy£(x,y)G(y,x') = <5(x-x') 



(29) 



The one-loop contribution to the free-energy density, rep- 
resenting the change in free energy due to fluctuations of 
the salt ions, can then be formally written as 



(3AF 



1 



log det K — \ log det 



-v; 



(30) 



To evaluate f3AF explicitly, we will employ the following 
mathematical identity for a general operator X: 



S log det X 



TvX^SX 



(31) 



We will first differentiate the free-energy density in Eq. 
(|30|l with respect to Ib- Using Eq. I|31|) . we can express 
the functional determinant of the kernel K(x, y) in terms 
of its operator inverse, the Green's function G(x, y). In 
the process, we need to make provisions to subtract out 
the (infinite) self-energy in the resulting expressions. In 
the form written in Eq. I|29[) . we can regard the Green's 
function G(x, y) as a fluctuating potential at point x 
generated by a test charge around which there exist fluc- 
tuating charges. Thus, the Green's function contains in- 
formation about correlations of the system, and in par- 
ticular, can be interpreted as the electrostatic interaction 
between two test charges located at x and y in the pres- 
ence of fluctuating ions in the bulk. We will now present 
the details on the computation of the Green's function. 



Construction of the Green's function 



where we have made a shift of variable u — z + zq and 
labelled the homogeneous solution h{u) and defined the 

quantity a 2 — \ (l + . The two independent homo- 
geneous solutions are: 



e -2a K u , 1 



hi(u) 
h 2 {u) = e 2aKU 



coth ku 
f 2a 
coth ku 
2a 



(35) 
(36) 



To construct the Green's function G(z, z'), we split space 
into three distinct regions and define in each of the region 
(Fig.EJ: 



G 



< 



1 > 



z' 



FIG. 3: Definition of the Green's function in three different 
regions: z < oo, < z < z' , and z > z . See Eqs. H37H . 138L 



To determine the Green's function, we use our action 
in Eq. (fTUf) to compute the kernel K(x, y), and we can 
write the representation of Eq. I|29l) in real space: 

{ - V 2 + K 2 cosh [<p (z)] } G(x, x') = Ib £(x — x') . (32) 

We note that the Green's function G(x, x') contains lat- 
eral invariance and can thus be Fourier-transformed in 
the directions perpendicular to the surface. Upon sub- 
stituting the saddle-point solution (po(z) [Eq. (j20J)] into 
the kernel, we obtain 



du 2 



+ q 2 + k 2 + 2k 2 csch 2 K (z + z Q )) 



G(z,z';q) 



= £ B 6(z-z'), 



(33) 

where q 2 = q 2 + q 2 . In order to solve Eq. (pTTTf) . we have 
to obtain the solution to the corresponding homogeneous 
equation: 



d 2 



4a 2 + 2 csch 2 , 



h(u) = , (34) 



G>(z,z') = A(z')h 1 {z) 1 z>z'(37) 
G<(z,z') = 5(z')Mz) + G(z> 2 (z) , 

< z < z' (38) 

G-(z,z') = D(z')h 1 (-z) , -oo<z<0 (39) 

To determine the coefficients A(z'), B(z'), C(z'), and 
D(z'), we impose the following boundary conditions: 



G_(0,z') = G<(0,z') 



<9G_ 



(40) 

dz 12=0 = ~dz~^ (41) 
G<(z,z')U=z' = G>(2,z')U=*' (42) 
dG< . 5G> . 



dz 



Or: 



(43) 



The first three conditions in Eqs. (gUJ), (gTJ, and j4^|) 

simply reinforce the continuity of the Green's function 
G(z, z') across the boundary z = 0. The presence of a 
charged-surface boundary at z = has already been in- 
corporated in the form of the mean-field solution <£o{z) 
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[Eq.lJ^njl], which we used in obtaining Eq. There- 
fore, the Green's function, which represents the correla- 
tion between two charges at distances z and z' from the 
charged surface, along with its derivative should be well- 
behaved and not contain any singularity. The last condi- 
tion in Eq. (|43(l represents the location of a test charge 
at z' as encoded by the delta function in Eq. (|32(l . So far 



in specifying the boundary conditions in Eqs. I|40|) - i|43|) 
we have only considered the case when z' > 0; results 
for the case of z' < can be deduced from symmetry. 
After some algebraic computation, we arrive at the fol- 
lowing expressions for the Green's function G(z, z') for 
the region < z < z': 



G(z,z';q) = ^ ^ n 2 + q 2 



-2 aK (z'-z) L ! coth«(z' + zo) \ / coth k(z + z Q ) 



2a 



2a 



+ M{q)e- 2aK(z ' +z) (l I cothK ( z ' + z °) \ ( 1 | coth k{z + z ) 
! \ 2a 1 \ 2a 



where 



M(q) 



V 1 + M) 2 



v/(<zA) 2 + (kA) 2 + y/1 + (kA) 2 ) 1 + V / (('?A) 2 + (k\) 2 ) (1 + (kA) 2 ) + (nX) 2 + (qX) 2 



(44) 
(45) 

(46) 



The Green's function for the other regions (G> and 
G-) can also be computed in a straightforward man- 
ner. For our purposes, we are only interested in the 



Green's function in the case of z = z', which represents 
self-correlation. G(z, z; q) takes the following expression: 



G{z,z;q) - 



coth 2 k(z + zq) 
Aa~ 2 



M(q)e 



coth k[z + Zq ) 
2^ 



(47) 



To understand the expression of the correlation self- 
energy, we see that in the n — > limit, the A / J(q)-factor 
in Eq. (|47|l reduces to the expression corresponding to 
the case of a single c harged plate with only counterions 
present in the system [26J, |5jJ • One might expect that the 
k — ► limit should correspond to a system with no ions. 
However, the condition of overall neutrality in the sys- 
tem necessitates the presence of counter-ions even when 
there is no salt present in the system. 

The self interaction is defined as the partially Fourier- 
transformed Green's function in Eq. (|47|l : 



G s (z, z) 



d 2 q 
(2n) 2 



G{z,z;q) - y~ 



(48) 



where we have subtracted out the infinite Coulomb self- 
energy at the end. The plot of a rescaled G s (z) is shown 
in Fig. As apparent in the figure, the self-energy is 
negative and indicates a decrease in the local free energy 
of the particles due to the screening of salt ions. This 



extends the same conclusion previously found in systems 
containing only counterions 48] . 



Interaction between Two Test Charges above 
the Surface 



As mentioned previously, the Green's function G(x, y) 
can be interpreted as the electrostatic interaction be- 
tween two test charges located at x and y in the presence 
of fluctuating ions in the bulk. Using the expression of 
the Green's function in Eq. (|47|l . we can revert back to 
physical space in the x- and y-directions and consider 
the interaction of two point charges separated by a dis- 
tance r from each other, and both at a distance z from 
the charged plate. The real-space correlation self-energy 
relates to Eq. I|47|l via a Hankel transform: 



G(z,r) 



(49) 
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FIG. 4: Plot of the rescaled partially transformed Green's 
function -^j—G a (z, z) as a function of kz at x — kA = 1. 



where Jo(qr) is the Bessel function of the first kind of 
order zero. When the two charged particles are far away 
from the charged surface (z — > oo), the interaction can 
be shown to reduce analytically to the Yukawa form: 



G(r) - £e- 



(50) 



Thus as one would expect, at far enough distance the two 
point charges act as if though there were no charged plate 
at all. On the other hand, for the two point charges that 
are close to the charged surface (small z), we can numer- 
ically integrate Eq. H49|) and plot the form of resulting 
interaction for various values of z, as shown in Fig. J5J). 

In Fig. JSJ, we compare the interaction between two 
point charges located in close proximity to the charged 
plate with the usual Yukawa-type interaction experienced 
by the same charges in the absence of or far away from 
the plate. As the distance z from the charged plate de- 
creases, we observe that the corresponding repulsion be- 
tween the two point charges of the same sign is dramat- 
ically suppressed due to density fluctuation. Therefore, 
we conjecture that when such other factors as Van der 
Waals interaction are taken into account, this decrease in 
electrostatic repulsion may lead to overall qualitatively 
different behavior for the two point-charges, such as at- 
traction between like charges (charge inversion), with an 
appropriate range of system parameters and in the region 
close enough to the charged plate. 



Fluctuation Contributions to Free-Energy 
Density 



FIG. 5: Plot of the rescaled correlation self-energy 
-j^G(r, kz) as a function of the scaled distance kz from the 
charged plate. We use x — kA = 0.1. We see that density 
fluctuation of the salt ions reduces the repulsion experienced 
between two point charges near the plate. 



functional determinant of the kernel K: 



(3AF = - log det K - - log det 



-V 2 



(51) 



The term proportional to det — V 2 ; comes from the 

normalization factor in the partition function and sim- 
ply subtract out the divergence associated with the self- 
energy of charged particles; a problem well-known even in 
the context of classical electrostatics. Physically, analo- 
gous to the mean- field result, the one- loop free energy 
consists of a volume contribution /3Afv and an area 
contribution (3AJa. For the kernel K in Eq. we 
see that all terms except for the last are translation- 
invariant. Therefore, when integrated over all space, 
these translation-invariant terms would produce factors 
proportional to the overall volume of space and would ex- 
actly constitute the volume contribution to the one-loop 
free energy. The remaining contribution has a depen- 
dence on the distance z from the plate and thus consti- 
tute the area contribution to the one-loop free energy. 

Therefore, to compute the one-loop correction to the 
free energy due to the presence of a charged plate, we 
must first isolate this translation-invariant piece of the 
kernel from the rest that actually depends on the ions' 
distance from the charged plate. This translationally in- 
variant piece in the kernel K corresponds precisely to 
what one gets in the Debye-Hiickel approximation for the 
salt case, where the electrostatic potential ip is linearized: 



Armed with the Green's function in Eq. Q47|). we can 
now calculate the fluctuation-induced correction to the 
free-energy density, formally expressed in terms of the 



K 



DH 



= -v 



<5(x - x') 



(52) 



The corresponding Green's function for the Debye-Hiickel 
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kernel is simply the Yukawa potential: 



G DH (z,z';q) = 



-1olk\z—z' I 



2^ 



2 + q 2 



(53) 



which shows explicit translational invariance. To deter- 
mine this volume contribution to the one-loop free en- 
ergy, we can simply consider the change in the internal 
energy density when a particle of charge e is brought 
into an electrolyte solution; the modification of this self- 
energy is exactly given by half the difference between the 
Yukawa potential and Coulomb potential in the limit of 
zero separation |57j |. Using standard thermodynamic re- 
lation between the internal energy and the free energy 
and the condition that Afv — > as T — > oo, we readily 
obtain the following result: 



1 



log det Kkh — \ log det 



12tt 



-V . 



(54) 



This is the classical result obtained by Debye and 
Hiickel 6l|. Thus, subtracting off the volume contribu- 
tion, we have the following expression for (3A/a' 



i 



log det K log det Kr>H 



(55) 



To proceed, we differentiate Eq. (|3U|) with respect to 
l B and apply the identity in Eq. (|3*T)l . As a result, the 
interesting part of the one-loop correction to the free- 
energy density is 



= k{l WfL ^^('.'iff)]^^ 00 ^^)]-/ dz[G DH (z,z;q)]^-[ K 2 ]+...} , 

(56) 



where A has the dimension of area. At this point, we 
need to make one more adjustment to render the expres- 
sion finite: we need to subtract off by the part in the 
Green's function that corresponds to the free-space case 



where there are no ions present, Go(q) = £ B /2 \q\, from 
each of the Green's functions in Eq. (|56|l . Upon further 
simplification, we can express the one-loop correction to 
the free-energy density in terms of three integrals: 



df3Af A 2 8k f d 2 q trr _ , 4 k 2 dX f d 2 q ^ 

1 {l 1+ 2I 2 }-— - TTTTo 1 ^ (57) 



di B Ib di B J (2tt) 2 l ^ i B y/l + (k\) 2 dl B J (2n) 2 

where 

poo 

11 = K dz [G(z, z;q) - G DH {z, z;q)} , (58) 

Jo 

f°° 2 

12 = k / dz [G(z, z; q) — Gdh (z, z; q)] csch k (z + zq) [1 — k(z + s) coth«:(z + zq)] , (59) 

Jo 

/•OO 

13 = K I dz [G(z, z; q) — Gq(z, z; q)] csch 2 K(z + zo)cothn(z + Zq) (60) 

Jo 

I 

It remains now to evaluate these integrals. We will one-loop correction, induced by salt density fluctuation, 

present the details of which in the Appendix. For now, we to the free-energy density in terms of the dimensionlcss 

simply state that after massaging these integrals further parameter x — kX: 
and integrating over l B , we would finally arrive at the 
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- -gfl (61) 

where 

VTT^ 2 



Kix) = ^ + v / 3 + 2x 2 -a; 4 tan" 1 — ~ = + 2x ( ^1 + £ 2 - x] + x 2 log ■ 

2V3 2x + VI + -^ 2 ^ ; 

1-x 2 



log[l + 7 _). (62) 



We note that the fluctuation-induced correction to the 
free-energy density has the overall effect of lowering the 
total free-energy density. This is in line with our expecta- 
tion for performing perturbation about a stable ground 
state. It also makes sense that the one-loop correction 
would lower the free-energy density since mean-field the- 
ory usually overestimates repulsive interactions due to 
the neglect of correlations. A plot of K(x) is shown in 
Fig. . We see that the one- loop correction to the free 
energy density saturates quickly when x > 2 in the high- 
salt limit. Therefore, in agreement with previous finding, 
in the regime where the interaction between salt ions at 
their mean separation is larger than thermal energy, we 
can simply define a renormalized screening length in the 
mean-field Poisson-Boltzmann description that results in 
a renormalized free-energy density. And the difference 
of which from the bare mean-field free energy does not 
depend sensitively at all on the exact value of the screen- 
ing length in the high-salt regime. However, fluctuation 
effects become important in systems with low to interme- 
diate salt-density characterized by < x < 2, and whose 
correction to the free energy density depends sensitively 
on the salt density as well as their charge valency. 



7C(x) 




FIG. 6: Plot of lC(x) as a function of x = k\. 



V. CONCLUSION 

We have formulated in a field-theoretic language the 
problem of a charged plate in a generally asymmetric 
salt solution containing both counterions and coions. We 
have obtained the non-linear Poisson-Boltzmann equa- 
tion as a saddle point of our field-theory action. Spe- 
cializing to the case of symmetric salt ions and using the 
full nonlinear mean-field solution, we systematically in- 
corporate the effect of fluctuation in ion densities and 
computed the one-loop correction to the free-energy den- 
sity in closed form by integrating a Gaussian in the action 
and using the techniques of Green's function. We have 
found that for systems with low to moderate salt density, 
fluctuation corrections to the free-energy depends sensi- 
tively on the salt concentration as well as their charge 
valency. We have also considered the interaction of two 
point charges near the charged plate and found that den- 
sity fluctuation suppresses repulsion between two point- 
charges. Therefore fluctuation in charge densities reduces 
the screening effect of the salt ions as well. 

Consequently, we see that the mean-field Poisson- 
Boltzmann approach may have wider applicability than 
we originally anticipate. Even in systems with high 
salt density or high-valency ions, the incorporation of 
their density fluctuation actually offsets the supposed 
stronger correlation between these ions, for which one 
would naively expect to invalidate the mean-field ap- 
proach. Thus in principle, using a suitably chosen "renor- 
malized" screening length, mean-field predictions can be 
made to fit properties of these charged systems in the 
high salt limit. We note that in contrast to the pertur- 
bation carried out in the strong-coupling regime [48| , the 
one-loop corrections that we have obtained are propor- 
tional the the surface charge density. Therefore, a high 
charge density for the plate should create a similar effect 
seen the high-salt limit. 

So far we have only considered the effects of electrolyte 
fluctuation on the electrostatic interaction. In reality, 
other factors such as the van der Waals forces are also 
essential to the physics of these systems. Though often 
practiced in calculations, the separation of electrostatic 
interactions from van der Waals ones is rather artificial, 
since quantum-mechanically it is the fluctuation in the 
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local electrostatic potential that gives rise to these at- 
tractive dispersion forces. It has been shown [EH. l5ij l6(i| 
that the zero-frequency van der Waals interactions are al- 
tered significantly by the presence of salt ions, while the 
nonzero-frequency dispersion can be effectively taken as 
independent of ionic densities since they do not respond 
to high-frequency fields. In regimes where linear super- 
position is applicable, our methodology and results for 
the one-plate case can be useful in studying the static 
contribution to these attractive dispersion forces. Our 
expression for the one-loop free energy should also pro- 
vide a more accurate description of the electrostatics in 
the study of phase transitions of a variety of systems that 
involve charged components in aqueous solutions. 

In conclusion, it is our hope that our solution would 
provide some insights as to why the mean-field Poisson 
Boltzmann approach seems to work in situations where 
the assumption of mean field should not be applica- 
ble [ill Il6j . It is our further hope that by means of 
such techniques as the Derjaguin approximation and 
superposition principle, our solution can be applied to 
provide better approximation even to the more general 

I 



where 



To integrate z over the range [0,oo), we first substitute 
Eqs. (|17J) and O into Eqs. (|X2|) - (|X4|) . and we repeat- 
edly apply integration by parts to obtain the following 



cases involving spherical symmetry or multiple charged 
plates in salt solution. The system of two charged plates 
in salt solution is of particular interest given the possible 
existence of attractive electrostatic forces between them. 
Our results for this case will be published elsewhere [§2 • 
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APPENDIX A: EVALUATION OF INTEGRALS 

In the appendix, we will present a few details on the 
evaluation of the integrals Ii , I2 , ^3 found in section IV- 
B, as follows: 



(Al) 

(A2) 
(A3) 
(A4) 



results, in terms of x — k\ and a = 




d(3Af A 2l^!L f d2q OTI 4 ^ d\ f d 2 g 

< u ' ~ B ~dT B J (2tt) 2 1 1 + 21 ~ Ib Vl + («A)2 9^ J 



d£i 



POO 

= k dz [G(z,z;q) - G DH (z,z;q)] , 
Jo 

/>oo 

= k I dz [G(z, z; q) — Gdh(z, z; q)] csch 2 K (z + zq) [1 — n(z + s) coth«;(z + zq)] , 



k dz [G{z, z; q) — Go(z, z; q)] csch k(z + zo)cothK(z + zo) 
Jo 



I2 = — 



-A{q) 
Ib ( 



4 ^ 2 V 



2 a 



1+ iyTT^ + 1 

4a \ ax Aa 



-l)-A(q) 



1 / v 7 ! + x 2 



- 1 



1 1 



4 x\J\ + x 2 



Mi) 

4x 4 



B(q) 



1 VT+. 

Act a; 3 
I 



1 



16a 2 x 2 



(A5) 
(A6) 
(A7) 



where 



A(q) 



tB K 



2q 2 ^Jn 2 + q 2 



B(q) = ^V^+q*M(q). 



(A8) 
(A9) 



To compute the momentum integrations, we use the stan- 
dard method of converting to polar coordinates in Fourier 
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space. Repeated application of partial-fraction decompo- 
sition in integrals involving M(q) and doing some alge- 



J 



braic manipulations lead to the following expressions: 



(27T) 



d 2 q 

W) 2 



47r V \/l + a; 2 



8tt 1 + x 2 



dy l + y 



dy 2 + y 



8irx 2 16nx 2 X y/l + x 2 J ^ 1 + V j^i + y + y 2 
I 



(A10) 
(All) 



where y = -j£=g. Using Eqs. I|A10|) and ilAllfl . we can and second term of Eq. (|A1|) . Upon further simplifying, 

now integrate away £ b on both sides of Eq. (JXTJ and per- we obtain the following intermediate expression for the 

form the k- and A-integrations, respectively, in the first one-loop free-energy density: 

I 



/3AjU 



4ttA 2 



where 



\J \ + x 2 — x^j — sinh 1 a; 



2U 



x dx 
1 + x 2 



4ttA 2 



2nX 2 ~ VT +^ + h lB > 



dx 1 
x 3 1 + x 2 



dy l + y 

y j^i + y + y 2 ' 

dy 2 + y 
l + y + y + y 2 



(A12) 

(A13) 
(A14) 



Evaluation of the remaining two integrals in Eqs. I|A13(I to the final expression of the one-loop correction to the 
and l|A14|l are straightforward. Substituting the subse- free-energy density in Eqs. I|tjl|) an d (|62|l . 
quent results back in Eq. (| A 12(1 and simplifying lead us 
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